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|We bound the rate of convergence to stationarity for a signed generalization of the 
Bernoulli-Laplace diffusion model; this signed generalization is a Markov chain on the 
homogeneous space (Z 2 I S n )/(S r x S n - r ). Specifically, for r not too far from n/2, 
we determine that, to first order in n, -jnlogn steps are both necessary and sufficient 
for total variation distance to become small. Moreover, for r not too far from n/2, we 
show that our signed generalization also exhibits the "cutoff phenomenon." 

1. Introduction. Consider the classical Bernoulli-Laplace model for the diffusion 
of gases through a membrane, in which at each step two randomly chosen balls from different 
urns are switched. How many steps does it take for this process to achieve near-randomness? 
This question was answered by Diaconis and Shahshahani (1987). Suppose that the balls 
also have charges and that, at each step, the two balls are not only switched, but their 
charges are also possibly changed. How many steps does it take for this process to achieve 
near-randomness? This is the question that we consider. 

For the Bernoulli-Laplace model, Diaconis and Shahshahani (1987) obtained bounds on 
the rate of convergence to stationarity. Similarly, in preparation for our main result, we 
bound the rate of convergence to stationarity for a variant of the classical Bernoulli-Laplace 
diffusion model; this variant is also a Markov chain on the homog eneous space S n j \S r x S n — r ). 
Specifically, for r not too far from n/2, we determine that, to first order in n, \n\ogn steps 
are both necessary and sufficient for total variation distance to become small. 

We then bound the rate of convergence to stationarity for a signed generalization of our 
variant of the Bernoulli-Laplace diffusion model; this generalization is a Markov chain on 
the homogeneous space (Z 2 I S n )/(S r x S n - r ). Specifically, for r not too far from n/2, we 
determine that, to first order in n, ^nlogn steps are both necessary and sufficient for total 
variation distance to become small. Moreover, for r not too far from n/2, we show that our 
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signed model also exhibits the "cutoff phenomenon." We also examine a slight variant of this 
signed model, establishing upper and lower bounds on its rate of convergence to stationarity. 

2. The Bernoulli— Laplace Diffusion Model. 

2.1. Introduction. We now review the Bernoulli-Laplace model for the diffusion of 
gases through a membrane. This is done in preparation for Section ||| where we extend the 
results of this section to a signed generalization of the Bernoulli-Laplace model. Imagine two 
racks, the first with positions labeled 1 through r and the second with positions labeled r + 1 
through n. Without loss of generality, we assume that 1 < r < n/2. Suppose that there are 
n balls, labeled 1 through n, each initially placed at its corresponding position. 

At each step, a ball is chosen at random from each rack and the two balls are switched. 
Then, if necessary, the balls on each of the two racks are permuted so that their labels are 
in increasing order. The state of the system is completely determined, at each step, by the 
labels of balls on the first rack. This switching procedure thus induces a Markov chain on 
the set of all (™) subsets of size r from a set with n elements. 

Let K be the subgroup of S n which permutes the first r indices among themselves and 
the last n — r indices among themselves, but does not commingle these two sets of indices. 
Thus K = S r x S n - r . Notice that K is the stabilizer of the initial configuration of the 
process. The switching procedure thus also induces a Markov chain on the homogeneous 
space X = S n / (S r x S n - r ). The set X may thus be identified with the set of all (™) subsets 
of size r from a set with n elements. 

Let T be the set of all transpositions in S n . Also let T\ be the set of all transpositions in 
K and let T 2 be the set of all transpositions in S n \ K. Thus T = T\ U T 2 . At each step, the 
process defined above chooses a random element of T 2 . 

But notice that before the two balls to be switched have been chosen, the balls in the first 
rack may be permuted among themselves and the balls in the second rack may be permuted 
among themselves, without changing the state of the system following the switch. Thus, at 
each step, the procedure actually chooses a random element of T 2 K := {tk G S n : r G T 2 
and k G K}. Notice that each element of T 2 K can be uniquely written as tk, where r G T 2 
and K G K. 

This Bernoulli-Laplace diffusion model may be modeled formally by a probability measure 
P on the symmetric group S n . Since \T 2 K\ = \T 2 \ ■ \K\ = r(n — r) ■ r\(n — r)\, we may thus 
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define the following probability measure on the set of all permutations of S n : 

P(tk) := — -— - -t where tk G T 2 K and 

r\n — r)r\[n — ry. (2 11) 

P(7r) := otherwise. 

Since there are n! elements in S n , the uniform probability measure on the set of all permu- 
tations of S n is given by 

U{ir) :=^- for every tt G S n . (2.1.2) 
n\ 

The following result, which is Theorem 2 of Diaconis and Shahshahani (1987), establishes 
an upper bound on both the total variation distance and the £ 2 distance between P* k and 
U, where P* k is the probability measure on the homogeneous space X = S n /(S r x S n - r ), 
induced by the convolution P* k of P with itself k times, and U is the uniform probability 
measure on X. (Homogeneous spaces and probability measures induced on them are discussed 
in Section |2.2| .) It should be noted that \\P* k — U\\tv is the total variation distance to 
uniformity after k steps for the Markov chain on X defined above, started at the chosen 
initial configuration. 

Theorem 2.1.3. Let P and U be the probability measures on the symmetric group S n 
defined in $2. 1 . 1\ ) and (\2.1.3j ), respectively. Let P* k be the probability measure on the homo- 
geneous space X = S n /(S r x 5* n _ r ) induced by P* k and let U be the uniform probability 
measure on X . Let k = ~r (l — -J (logn + c). Then there exists a universal constant a > 
such that 

||P**-£/||tv < |C) V2 ||^-^ll2 < ae~ c/2 for all c> 0. 
In the following sections we present the results that were needed to prove this theorem 



and which are used to prove analogous results in Section p.5| and in Section y. In Section [L2 



we present the definitions and results necessary to study Markov chains on homogeneous 



spaces. In Section gT3| we decompose the representation L(X) (defined in Section |272| ) of S n 



into its irreducible constituents. In Section 2.4 we calculate the Fourier transform of the 



probability measure defined in ( 2.1.1]) , using a procedure different from that of Diaconis and 



Shahshahani (1987), and show how it was used to prove Theorem 2.1.3 . In Section 2J5 we 
perform a similar analysis on a variant of the classical model; this variant will be generalized 
in Section El. 
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2.2. Markov Chains on Homogeneous Spaces. We now present basic properties 
and results regarding homogeneous spaces and Markov chains on them. A more detailed 
introduction to this subject may be found in Chapter 3 of Diaconis (1988). 

An action of a group G on a set A is a mapping from the Cartesian product G x A to A, 
with the image of (g, x) under this action being denoted by gx, which satisfies (i) ex = x for 
the identity element e G G and all x G A, and (ii) (g2gi)x = g2{g\x) for all g\,g2 G G and 
x G A. For an element x G A, the set of elements {g G G : gx = x} is called the stabilizer 
of x\ the stabilizer of an element x G A is a subgroup of G. A group G acts transitively on 
a set A if, for every Xi,X2 G A, there exists some g G G such that gxi = X2- A set with a 
group acting transitively on it is called a homogeneous space. 

Suppose that G acts transitively on a set A and that A is the stabilizer of some fixed 
element xq G A. The group G acts on the left coset space G/K by sending (g, hK) to {gh)K 
for all g G G and /i A G G/K. The set A and the left coset space G/K are isomorphic under 
this action. We may thus identify A with a set of left coset representatives {x , x%, . . . , x m _i} 
of K in G, where xq = e G G and ;Ej G G for 1 < i < m — 1. 

A probability measure P defined on a finite group G induces a probability measure P on 
the set A = G/A by defining P{xj) := P(xiK) for < z < m — 1, where a^A := {xj/c G G : 

G A}. The transition matrix P of the Markov chain induced by the probability measure 
P is defined, for the element at the intersection of the row corresponding to Xi G A and the 
column corresponding to Xj G A, by 

Pxi.Zj := P^XjKXi ) 

where XjKx^ 1 := {xjkx^ 1 G G : G A}. According to Lemma 3 in Section F of Chapter 3 
of Diaconis (1988), P is doubly stochastic. Furthermore, the fc-step transition probabilities 

~k 

P for the Markov chain satisfy 



and therefore, since x := e, 



■k 



= P* k (xjK) = P*H 



This confirms the statement concerning \\P* k — U\\tv that immediately precedes Theo- 
rem 



2.1.3 



Let V be a finite-dimensional vector space over the complex numbers and let GL( V) be the 
general linear group of isomorphisms of V onto itself. Let p : G — > GL(V) be an irreducible 
representation of G with character \p an d dimension d p . Suppose that H is a subgroup 
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of G and that p\ and pi are representations of G and H, respectively. Let p\ be the 
representation of H known as the restricted representation and let p2 be representation 
of G known as the induced representation. 

Let L(X) be the set of all (complex- valued) functionals defined on X. The group G acts on 
L(X) via the definition (gf)(x) := f(g~ 1 x) for g G G, x G X, and / G L(X). Since for fixed 
g G G this action is a bijective linear mapping of L(X) onto itself, L(X) can also be regarded 
as a representation of G. In fact, as a representation, L(X) is isomorphic to po the trivial 
representation p of .fT induced up to a representation of G. According to Theorem 2 in 
Section 1.4 of Serre (1977), every representation of a group G is the direct sum of irreducible 
representations of G. Thus L(X) is the direct sum of irreducible representations of G. 

According to Theorem 6 in Section 2.5 of Serre (1977), the characters xi, X2, • • • , Xs of the 
irreducible representations of a group G form an orthonormal basis for the Hilbert space of 
class functions on G with respect to the inner product defined by 

(i ) ,x) G TQ^Mxig)- 

1 1 gaG 

For any irreducible representation p of G, define 



m p := (x p ,Xl(x)) g - 

Thus m p is the multiplicity of the irreducible representation p of G in the decomposition of 
L(X) into its irreducible constituents. 

The following useful result relating induced and restricted characters is the Frobenius 
reciprocity formula, which is Theorem 13 in Section 7.2 of Serre (1977). 

Lemma 2.2.1. Suppose that H is a subgroup of G and that if> and x ore characters of H 
and G, respectively. Then 

(^1h,x) g = (*P,XI%) H , 
where the inner product on the left (resp., right) is calculated in G (resp., H). 



It thus follows from Lemma 2.2.1 that 



™ P = (x p ,Xl(x)) g = (Xp,Xpo1k) g = (XpIk,Xpo) k - 

Thus m p is also the multiplicity of the trivial representation po of K in the decomposition of 
the restriction of p to K. Furthermore, we have the following, which is Lemma 5 in Section 
F of Chapter 3 of Diaconis (1988). 
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Lemma 2.2.2. The multiplicity m p of an irreducible representation p : G — > GL(V) of 
G in the decomposition of L(X) into its irreducible constituents equals the dimension of the 
space of K -fixed vectors in V , i.e., 

m p = dim{t> G V : p{k)v = v for all k G K}. 

For notational purposes, define I(p) := I © 0, where I is the m p -dimensional identity 
matrix and is the (d p — m p ) x (d p — m p ) zero matrix. Thus I(p) is the d p x d p matrix 

g g . Notice that tv(l(p)I(p)*\ = m p . The preceding lemma leads to the following result, 
which will be useful in the calculation of the Fourier transform in Sections 12.41 and 



Lemma 2.2.3. Suppose that a finite group G acts transitively on a finite set X and that 
K is the stabilizer of some fixed element Xq G X. Also suppose that p : G — ► GL(V) is an 
irreducible representation of G. Then there is an orthonormal basis ofV such that 

= fa)- 



Proof. Lemma |2.2.2j asserts that m p is the dimension of the space of i^-fixed vectors 



in V. Choose an orthonormal basis in V such that the first m p basis vectors are i^-fixed. It 
then follows from Theorem 1 in Section 1.3 of Serre (1977) that, under this choice of basis, 
the representation p splits as V = V\ © V 2 , where dim V\ = m p and dim V 2 = d p — m p . 
Thus p ljf= pi © p 2 , where p 1; as a representation on V±, is the direct sum of m p trivial 
representations on K and p 2 is a representation on V 2 . 

Notice that — - p(k) is the Fourier transform of the uniform distribution on K at the 

' ' fceif 

representation p of K. Furthermore, for any finite group, the Fourier transform of any 
probability distribution at the trivial representation is 1 and the Fourier transform of the 
uniform distribution at any nontrivial representation is the zero matrix. It is from this that 
the desired result follows. □ 

When m p < 1 for every irreducible representation p of G, the decomposition of L(X) is 
called multiplicity free. In such a case the pair (G, K) is called a G elf and pair. Thus if (G, K) 
is a Gelfand pair and p is an irreducible representation of G occurring in L(X), then I(p) is 
the dp x dp matrix whose only nonzero entry is a "1" in the (1, 1) position. A more detailed 
introduction to Gelfand pairs may be found in Diaconis (1988) or Macdonald (1995). 

A probability measure P defined on a finite group G is called right K-invariant if P{gk) = 
P{g) for all g G G and k G K. There is an analogous definition for left i^-invariance. A 
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probability measure P is called bi-K-invariant if it is both left and right i^-invariant, i.e., 
if P(kigk,2) = P{g) for all g G G and all k\,k2 G K. Our probability measure P defined 
in (|2.1.1|) is clearly right .^-invariant. Notice for any tk G T2K and any K\ G K that 
kitk = (kitk^ 1 )(kik) with K\TK^y G T 2 and K\K G K. It follows that P is also left K- 
invariant and hence bi-i^-invariant. 

According to Theorem 9 in Section F of Diaconis (1988), if (G,K) is a Gelfand pair and 
P is a bi-i^-invariant probability measure defined on G, then the Fourier transform P(p) at 
any irreducible representation p of G is a constant times I(p), where I(p) is the d p x d p 
matrix whose only nonzero entry is a "1" in the (1, 1) position. 

We are now able to bound the distance to uniformity of a probability measure P induced 
on a homogeneous space X = G/K from a right i^-invariant probability measure P defined 
on a finite group G in terms of the Fourier transform of P by use of the following, which is 
the Upper Bound Lemma in Section F of Chapter 3 in Diaconis (1988). 

Lemma 2.2.4. Suppose that a finite group G acts transitively on a finite set X and that 
K is the stabilizer of some fixed element xq G X . Also suppose that P is a right K -invariant 
probability measure defined on G, that P is the induced probability measure on the set X , and 
that U is the uniform probability measure on X . Then 

\\P-U\$v < l\X\'\\P-U\\l = |]Trf p tr (P(p)P(p)*) 

p 

where the sum is taken over all nontrivial irreducible representations of G that occur at least 
once in L(X). 

Notice in the special case K = {e} that X = G / K = G and that the preceding lemma reduces 
to the Upper Bound Lemma for groups (see, e.g., Diaconis (1988), Chapter 3, Section B). 

2.3. Irreducible Representations in L(X). As a representation of S n , L(X) is isomor- 
phic to p fs" x Sn _ r , where p is the trivial representation of S r x S n - r . The decomposition 
of L(X) into its irreducible constituents is given by the following, which is Lemma 2.2.19 of 
James and Kerber (1981). 

Lemma 2.3.1. Let L(X) be the set of all (complex-valued) functionals defined on X = 
S n /{S r x S n - r ). Then, as a representation of S n , 



L(X) = p [n] ® P[n-l,l] © ••• © P[n-r,r] 
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where p\ n -i,i] is the irreducible representation of S n corresponding to the partition [n — i,i] of 
n, forO <i <r. Furthermore, the dimension of p\ n -i t i] is given by 




with [n — 0, 0] = [n] and the the usual conventions (^) = 1 and ( * ) = 0. 



Notice that the decomposition of L(X) is multiplicity free. Thus (S n , S r x S n - r ) is 
a Gelfand pair. So for any irreducible representation p occurring in the decomposition of 



2.4. Analysis of the Classical Model. In order to apply the Upper Bound Lemma 
( p.2.4| ) to our bi- .^-invariant probability measure P, we must now calculate the Fourier 
transform at each nontrivial irreducible representation of S n occurring in the decomposition 
of L(X). Our calculations are done with the aid of the following two lemmas. 

In the special case when P is a class function, it is a consequence of Schur's Lemma that 
the Fourier transform may be calculated easily by use of the following, which is Lemma 5 of 
Diaconis and Shahshahani (1981). 

Lemma 2.4.1. Suppose that p is an irreducible representation of a finite group G with 
character \ and that P is a class function. For each conjugacy class i, let Pi be the constant 
value of P on the class, let be the cardinality of the class, and let Xi be the constant value 
of x on the class. Then the Fourier transform of P is given by 



where d p is the dimension of p, I is the d p - dimensional identity matrix, and the sum is taken 
over distinct conjugacy classes. 

The following formulas, found in Section D of Chapter 3 and Section B of Chapter 7, 
respectively, of Diaconis (1988) are used to calculate the numerical value of the Fourier 
transform. 

Lemma 2.4.2. Suppose that p is an irreducible representation of S n corresponding to the 
partition [A] = [Ai, . . . , Xk] of n. Let r(A) := X[\]( T )/d[x] with r e S n . Then 



L(X), I(p) is the d p x d p matrix whose only nonzero entry is a "1" in the (1, 1) position. 



•p i=i 



r(A) 




with 1/m! := if m < 0. 
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We now calculate the Fourier transform at each nontrivial irreducible representation of S n 
occurring in the decomposition of L(X). Diaconis and Shahshahani (1987) did this with the 
aid of spherical functions. Our technique is different and was used by Greenhalgh (1989) and 
Scarabotti (1997) in their analyses of variants of the Bernoulli-Laplace diffusion model. 



Lemma 2.4.3. Let P be the probability measure on S n defined in (\2.1.J\) . Let p = p[ n ~i,i] 
(for some 1 < % < r) be a nontrivial irreducible representation of S n occurring in the decom- 
position of L(X). Then, in a certain basis, the Fourier transform is 

i(n — i + 



P(P) 



1 - 



I(P) 



r(n — r) 

where I(p) is the d p x d p matrix whose only nonzero entry is a " 1" in the (1, 1) position. 

Proof. Recall that pr n i is the trivial representation of S n . Thus we must calculate the 
Fourier transform for the other r nontrivial irreducible representations of S n occurring in the 
decomposition of L(X), which were found in Lemma |2.3.1| . Notice that, in the notation of 
Section pTTl , 

P{tk) P {tk) ' 



P(P) 



tk£T 2 K 



r(n — r) 



r{n — r)r\{n — r) 



r!(n — r)! 



Since T is a conjugacy class in S n , it follows from Lemma |2.4.1| that 

^p(r) = \T\r(p)I = \n{n - 1) r(p) I, 



where r(p) := Xp( T )/d p , with r G T C S n , and / is the <i p -dimensional identity matrix. 

Since L(X) = p$ |^ n , where po is the trivial representation of K, it follows that for any 
p occurring in the decomposition of L(X), p J.^ 1 is the direct sum of d p copies of po- Thus 
since T\ <Z K, we have 



[|n(n — 1) — r(n — r)\ I. 



£'00 - 

Now choose an orthonormal basis in V such that the first m p basis vectors are K- invariant, 



as described in Section [2.2| . It then follows from Lemma [2.2.3| that, in this basis, 

1 



r\(n — r)\ 



£p(«) 



Hp)- 
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Therefore, it follows from these results that 

Hp) - 



r \n{n — l)r(p) — \n{n — 1) + r{n — r) 



Hp)- 



r{n — r) 

We now calculate r(p) by use of Lemma [2.4.2 . For the ith nontrivial irreducible repre- 
sentation p = P[ n -i,i] of S n occurring in the decomposition of L(X), it follows that, for 
1 < i < r, 

1 



r(p) 



n(n — 1) 
1 



[{n-if - (n-i) + i 2 - 3i] 
[(n-i)(n-z-l) + i(i-3)] 



and hence 



n(n — 1) 

|(n — i)(n — i — 1) + — 3) — — 1) + r(n — r) 



1 



i(n — i + 1) 
r(n — r) 



r(n — r) 

T(p)- □ 



Hp) 



We have now established the results necessary to prove Theorem [2.1.3| . Recall from Sec- 
tion |2.2| that the probability measure P defined in Q2.1.1]) is bi-i^-invariant. By applying the 
result from Lemma [2.4.3| to the Upper Bound Lemma fl2.2.4p , we find that 

n? a -c'iiTv<i(;')ii^ i: -f , ii3 



P m p 



i(n + 
r(n — r) 



2k 



where the sum is taken over all nontrivial irreducible representations p = p\ n -i,i] of S n 
occurring in the decomposition of L(X). 



Since, for each of these irreducible representations, d p was determined in Lemma [2.3. 1| and 

p** -u\\% v <\{ n ^\\W k -u\\\ 



vn p = 1, it then follows that 



< 



iE 



i=i 



iE 



i=l 



n 
i - 1 

n 
i - 1 



exp 



r[n — r 



2ki(n + 
r(n — r) 



2k 



Thus, when k = ~r (l — -) (logn + c), 
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|p**-^lTv<!Qll^* fc -^ll 

r 



i=l 



n 
i - 1 



n ~ i{n—i+l)/n —d(n— i+l)/n 



A detailed discussion in Section 3 of Diaconis and Shahshahani (1987) determines the 
existence of a universal constant a > such that, for c > 0, 



IE 

i=i 



n 
i - 1 



n -i(n-i+l)/n e -ci(n-i+l)/n < ^ ^ 



(2.4.4) 



This completes the proof of Theorem 2.1.3| . 

Theorem p. 1.3 shows that k = \r (l — ^) (logri + c) steps are sufficient for the (normal- 
ized) £ 2 distance, and hence also the total variation distance, to become small. Diaconis and 
Shahshahani (1987) established a matching lower bound in the special case r = n/2. 

2.5. Analysis of a Variant of the Classical Model. In preparation for our analysis in 
Section |3], we now introduce a variant of the Bernoulli-Laplace diffusion model. Suppose that 
the balls and racks are as described in Section |2.1| . At each step of our variant of the classical 
model, two positions p and q are chosen independently and uniformly from {1,2, ... ,n}. If 
p 7^ q, switch the balls in positions p and q. If p = q (which occurs with probability 1/n), 
leave the balls in their current positions. Then, if necessary, permute the balls on each of the 
two racks so that their labels are in increasing order, as in the classical model. 

This switching procedure is exactly that of the random walk on the set of all possible 
positionings of n cards introduced in Section 1 of Diaconis and Shahshahani (1981). Fur- 
thermore, if 1 < p, q < r or (r + 1) < p, q < n, then the resulting state of the system is 
unchanged; this occurs with probability 1 — 2r ( T1 ~ r ) . This procedure induces a Markov chain 
on the same state space as the classical model, but slows down the process by a factor of 

one 



by not forcing two balls to be switched at each step. In light of Theorem p. 1.3 



2r(n— r) 

would certainly then expect |n(logn + c) steps to suffice for near-stationarity. We establish 
this conjecture as Theorem 2.5.3| . 

Notice that K = S r x S n - T is also the stabilizer of the initial configuration of this process. 



Let T, T\ and T2 be defined as in Section |2.1| . As with the classical model, before the two 
positions p and q have been chosen, the balls in the first rack may be permuted among 
themselves and the balls in the second rack may be permuted among themselves, without 
changing the state of the system following the switch. Thus, at each step, the procedure 
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actually chooses a random element of T 2 K := {tk G S n : r G T 2 and /t G X} with 
probability 2r( - n ~'' ) an d chooses a random element of if with probability 1 — 2r ( n ~ r > . 



As in Section [2.1| , this procedure may be modeled formally by a probability measure Q on 
the symmetric group S n . Since Q(tk) = 2r ^ 2 ~ r ^ P(tk) for all tk G where P is defined 

in (|2.1.1j) , we thus arrive at the following probability measure on the set of all permutations 
of S n : 



Q{k) :-- 

Q{tk) :-- 
Q(tt) :- 



n — 2r{n — r) 
n 2 r\(n — r)! 
2 



n 2 r\(n — r)! 
otherwise. 



where n G K, 
where tk, G T 2 K, and 



(2.5.1) 



For this probability measure Q, we now calculate the Fourier transform at each nontrivial 
irreducible representation of S n occurring in the decomposition of L(X), as was done in 
Section ITM for P. 



Lemma 2.5.2. Let Q be the probability measure on S n defined in ( \2.5.1\ ). Let p = p\ n -i t i] 
(for some \ <i <r) be a nontrivial irreducible representation of S n occurring in the decom- 
position of L(X). Then, in a certain basis, the Fourier transform is 



Q{p) 



2i(n 



1) 



n- 



Hp) 



where I(p) is the d p x d p matrix whose only nonzero entry is a " 1" in the (1, 1) position. 



Proof. Let P be the probability measure on S n defined in ( |2.1.1[ ). Notice that 



Q 



2r(n — r) 



rr 



U 



2r(n — r) 



K 



n 



where Uk is the probability measure on S n defined by Uk(k) :— t^t for all k G K and 
Uk{k) := otherwise. It then follows from Lemmas |2.2.3| and |2.4.3| , using the same basis, 
that 
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Q(p) = E QWM 



TT£Sn 



2r(n — r)\ , . . . 2r(n — r) , . , 

i - 1 = ; I E m^OpW + n2 E 



rr 



7T 



i - 2r(n r r) ^ jo) + 2r(n r r) p( P ) 



2r(n — r)\ 2r(rr — r) / z(n — i + 1) 

r[n — r) 



rr 



rr 



Hp) 



2i(n -i + 1) 



rr 



Hp)- □ 



The following result establishes an upper bound on both the total variation distance and 
the £ 2 dist ance between Q* k and U, where Q* k is the probability measure on the homogeneous 
space X = S n J (S r x S n - r ), induced by the convolution Q* k of Q with itself k times, and U 
is the uniform probability measure on X. 



Theorem 2.5.3. Let Q and U be the probability measures on the symmetric group S n 
defined in (\2. 5.1 ) and ( \2.1.2\ ), respectively. Let Q* k be the probability measure on the homo- 
geneous space X = S n /(S r x S n - r ) induced by Q* k and let U be the uniform probability 
measure on X . Let k = ^n(logn + c). Then there exists a universal constant a > such that 

\\Q*-U\\tv < i(") 1/2 ||Q^-?7|| 2 < ae~ c / 2 for all c> 0. 



This Markov chain is twice as fast to converge as the random walk on the symmetric group 
S n introduced in Section 1 of Diaconis and Shahshahani (1981). 



Proof. Notice that Q is clearly right 7^-invariant. (In fact, it is bi-TT-invariant.) By 
applying the result of Lemma 2.5.2 to the Upper Bound Lemma ( 2.2.4Q , we find that 



\\Q 



*k _ JJW 2 



TV 



< lO\\Q* k -u\\i = iE rf 



v p m p 



1 - 



2i(n-i + l) 



rr 



2k 



where the sum is taken over all nontrivial irreducible representations p = P[ n -i,i] of S n 
occurring in the decomposition of L(X). 
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Since, for each of these irreducible representations, d p was determined in Lemma [2.3. 1| and 



m p = 1, it then follows that 

\\^-u\\^<\Q)\\Q*-u\\ 

r 

= 1E 



1=1 



i=l L 



n 
i - 1 

n 
i-1 



2i(n 



2k 



(2.5.4) 



exp 



4fcz(n -i + 1) 



Thus, when k = jn{\ogn + c), 

IIQ 5 -c'llTv<i(;)IIQ 3 -c'l 



^iE 

i=l 

It then follows from ( |2.4.4| ) that, for c> 

r 

IE 



i - 1 



n -i( n — *+l)/ n g — ci(n— i+l)/n 



i=l 



z - 1 



^-i(n-i+l)/n e -ci(n-i+l)/n 



from which the desired result follows. □ 



Theorem |2.5.3j shows that k = \n (log n + c) steps are sufficient for the (normalized) 



distance, and hence the total variation distance, to become small. A lower bound in the 
(normalized) £ 2 metric can also be derived by examining (n — 1) (l — f) 2fc , which is the 
dominant contribution to the summation ( 2.5.4T ) from the proof of Theorem |2.5.3| . This term 
corresponds to the choice i — 1. Notice that k — \n (logn — c) steps are necessary for just 
this term to become small. 

A number of steps necessary for the total variation distance to become small is given by 
the following result. An analogous result was proved for the classical model (|2.1.1| ), in the 
special case r = n/2, in Theorem 1 of Diaconis and Shahshahani (1987). 

Theorem 2.5.5. Let Q and U be the probability measures on the symmetric group S n 
defined in \2. 5. 1\ ) and $2.1.2j , respectively. Let Q* k be the probability measure on the homo- 
geneous space X = S n / (S r x 5* n _ r ) induced by Q* k and let U be the uniform probability 
measure on X . Let n > 3 and 

k = min {ln(\ogn - c), (n - 2) log (^) + \{n - 2) log [l + \{n - 2)(l - {^f) e~ c ] } 
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be a nonnegative integer, with c > (and where log (^) is interpreted as +00 when r = n/2). 
Then there exists a universal constant a > such that 



iQ 1/2 ||Q**-[7|| 2 > ||Q**=-[7|| TV > 1- 



ae 



Proof. Recall from Lemma |2.3.1| that for L(X), the set of all (complex- valued) func- 
tionals defined on X = S n /(S r x S n - r ), we have the multiplicity-free decomposition 

L(X) = p [n ] © P[n-l,l] © ■•• © P[n-r,r], 



where P[ n -i,i] '■ S n — > GL (Vi) (say). So it follows from Lemma |2.2.2j that each Vi, for < % < 
r, has a unique nontrivial one-dimensional space of K- fixed vectors, where K := S r x S n _ r . 

The spherical function Sj, for < i < r, is defined to be the unique left i^-invariant function 
(i.e., vector) in Vi normalized so that Sj(xo) = 1, where xo is the left coset representative 
corresponding to e 6 S n . It follows from Exercise 17 in Section F of Chapter 3 of Diaconis 
(1988) that the spherical functions satisfy 

Si{x) = Tj-r\^2xi{xK) forxeX, 

where Xi is the character of p\ n -i,i], with [n — 0, 0] = [n]. 
Under the uniform measure U on X, notice that 

E u( s i) = ]x\J2 s i( x ) = Tx\^2 \K\^2xi(xK) = — { ^Xi^) = (Xi,Xo)s n , 
' ' xex ' ' xex ' ' neK ' 7reS„ 

where Xo is the character of the trivial representation pr n i. Thus since the irreducible char- 
acters of a group G form an orthonormal basis for the Hilbert space of class functions on 67, 
it then follows that E^j(so) = 1 and Ejj{si) = for 1 < % < r. 
Under the /c-fold convolution measure Q* k on X, notice that 



xex x&x ' ' kgk 



n—i,i\ 1 1 



where we use the fact that Q* h (x) = \K\ Q* k {xn) for each k e K. In particular, it follows 
from Lemma |2.5.2| that 

E^isi) = {I- If and E~ k (s 2 ) = (l - 
Define f(x) := \/n — 1 Si(x) for x & X. Then 
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E~ k (f) = v£=T(l- If and EfjU) 



0. 



In order to determine Var^(/) and Var~ fc (/), we must calculate E^(f 2 ) and E~ k (f 2 ). This 
is done with the aid of the following identity (which can be derived from formulas for s± and 
S2, as suggested in Diaconis and Shahshahani (1987)): 

1 



(n — 2r) 2 n 2 (r — l)(n — r — 1) 

n — 1 r(n — r)(n — 2) r(n — r)(n — l)(n — 2) 



Thus 



Var & (/) = (n - l)i%(s 2 ; 



n-_l _l (n-l)(n-2r) 2 F _ ^ ^ _,_ (n-l)n 2 (r-l)(n-r-l) 

r(n-r)(n-l)(n-2) ^C/ 



n-1 V(n-r)(n-2) -^C/( S l) + 



and 



Var~ fc (/) = (n-l)£~ fe (s 2 ) - (n - 1) (^(ai))' 

(n-l)(n-2r) 2 ^ 2^ 



n-1 
n-1 



+ 



r(n— r)(n— 2) 



+ 



(n-l)n 2 (r-l)(n-r-l) 
r(n— r)(n— l)(n— 2) 



2k 



1 , 4 (S5j)(n^) 2 , _ 2 y 
1 n 2 -(n-2r) 2 ^ nJ 



_ 3n-2 

n 2 — (n— 2r) 2 n— 2 



(1-1) 



2/,- 



Therefore, 



Var~ fe (/) 



1 



in 2 



(2.5.6) 



+ 



n-2 



n 2 — (n — 2r) s 



n — 2r 



n 



2 

1 - - 

n 



-k 



+ 



3n - 2 



(n - l)(n - 2) 



By elementary calculus, x < — log(l — x) < for < x < 1. Thus if k < \n (logn — c), 
with n > 3 and c > 0, then 



-2fc/(n-2) 



Q 



> (1 - I) 1/2 (I) 1/{n - 2) e ^ e -/in-2) > /J e c/2 
— \ n/ \nJ — V 27 



1 /2 

where we note that, for n > 3, (1 — i) (-) 



n l/2 / n l/(n-2) . . . c/(n-2W 1 

' ' ' is increasing and e / ( > > 1. 



Also notice that the first and third terms on the right in ( 2.5.6 ) are bounded by 
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3n - 2 



2 — T e C an d 



< 



3n 



(n-l)(n-2)" (| n ) (|n) 
when n > 3 and < c < logn. In order for the second term on the right in ( 2.5.6[ ) to be 
bounded above by e~ c we must have 



n — 2r 



n 



n 



< 



(n - 2)[n 2 - (n - 2r) 2 ] 
1 + 4^ 



(2.5.7) 



Since, when n > 3, 



n 



< 2fc/(n-2) 



it is sufficient for ( |2.5.7|) to have 



2A; 
n - 2 
i.e., to have 



< 21o£ 



n 



n — 2r 



fc < (n - 2) log 



n 



n — 2r 



log 



|(n-2)log 



2)[n 2 - (n- 2r) s 



1 + 



An 2 

(n - 2)[n 2 - (n - 2rf 
An 2 



In summary, if n > 3 and < c < log n, and if 

k = min |in(logn - c), (n - 2) log (^) + ±(n - 2) log [l + \{n - 2)(l - (^r) 2 ) e~ c ] }, 



then 



E^f) > \/4fe c/2 and 



Var~ fe (/) 



< 28e~ 



Now define A a := {x G X : |/(x)| < a}. It follows from Chebyshev's inequality that 
U(A a ) > 1 and that C?* fc (A*) < - r -^ = ^, provided < a < J^e c/2 . Then 



a 1 



fC) 1/2 HQ* fe -t/|| 2 > \\Q* k -u\\ 



vv 27 



TV 



^ 1 " — 

or 



28e~ 



| e</ 2 - a 



2 ' 



Choosing a = | y e C//2 shows that 
IC) 1/2 ||^-^||2 > HQ^-^Htv > 1 - 54e" c - 1512e- 2c > 1 - 1566 e - c , 
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which completes the proof. □ 

Theorem [2.5.5| gives a number of steps necessary for the total variation distance to become 
small in our variant of the classical Bernoulli-Laplace diffusion model. We now examine in 
rough terms how the value of r dictates the choice of k. For notational purposes, let 

f(n,r,c) := (n - 2) log (^) + ±(n - 2) log [l + J(n - 2) (l - (^) 2 ) e~ c ~ . 

Notice that f(n, r, c) is increasing in r. At one extreme, let r = n/2 (assuming for simplicity 
that n is even); then, since f(n,r,c) = +oo, fc is chosen to be |n(logn — c), matching the 
upper bound in Theorem [2.5.3| . At the other extreme, let r = 1; then it can be shown that 
^e~ c n < /(n, r, c) < |n (logn — c), and so k is chosen to be f(n, r, c) > ^e~ c n. In this case 
we find that order n (not order nlogn) steps are necessary, and indeed it is easy to show 
that order n steps are also sufficient. There is some value of r for which f(n,r,c) "crosses 
over" \n (logn — c). This occurs in the vicinity of 2r = n 1 / 2 . But even for 2r as small as n s , 



for any fixed 5 > 0, the value of k determined in Theorem [2.5.5| is of order nlogn 



3. A Signed Generalization of the Bernoulli— Laplace Diffusion 
Model. 



3.1. Introduction. We now extend (the variant in Section of) the Bernoulli-Laplace 
diffusion model to the case in which the balls also have charges (positive or negative). Imagine 
two racks, the first with positions labeled 1 through r and the second with positions labeled 
r + 1 through n. Without loss of generality, we assume that 1 < r < n/2. Suppose that 
there are n balls, labeled 1 through n, each initially placed at its corresponding position. 
Also suppose that each ball has a charge (positive or negative) and that initially each ball is 
positively charged. We refer to this as the signed Bernoulli-Laplace diffusion model. 

At each step, independently choose two positions p and q uniformly from {1,2, . . . , n}. 
If p ^ q, switch the balls in positions p and q. Then independently, with probability 1/2, 
change the charge of the ball moved to position p; and independently, also with probability 
1/2, change the charge of the ball moved to position q. Then, if necessary, permute the balls 
on each of the two racks so that their labels are in increasing order. If p = q (which occurs 
with probability 1/n), leave the balls in their current positions. Then, again independently 
with probability 1/2, change the charge of the ball in position p = q. 

We refer to the process described above as the independent flips model. A similar process, 
known as the paired flips model, is introduced in Section |3~j . 
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This switching procedure is exactly that of the random walk that was introduced in Section 
3.1 of Schoolfield (1999) in the special case of the hyperoctahedral group Z 2 I S n . The state of 
our signed Bernoulli-Laplace system is completely determined, at each step, by the ordered 
n-tuple of charges of the n balls 1,2, ... ,n and the labels of balls on the first rack. Our 
switching procedure thus induces a Markov chain on the set of all 2 n • (") ordered pairs of 
n-dimensional binary vectors and r-element subsets of a set with n elements. 

Let K be the subgroup of Z 2 I S n which permutes the first r indices among themselves and 
the last n — r indices among themselves, but does not commingle these two sets of indices. 
Thus K = S r x S n _ r . Notice that K is the stabilizer of the initial configuration of the 
process. The switching procedure described above thus also induces a Markov chain on the 
homogeneous space X = (Z 2 I S n )/(S r x S n - r ). The set X may thus be identified with 
the set of all 2 n ■ (™) ordered pairs of n-dimensional binary vectors and r-element subsets of 
a set with n elements. 

Let T be the set of all signed transpositions in Z 2 I S n . Also let 7\ be the set of all signed 
transpositions in K, let T 2 be the set of all signed transpositions in (Z 2 I K)\K, and let 
T 3 be the set of all signed transpositions in (Z 2 I S n ) \ (Z 2 I K). Thus T = T x U T 2 U T 3 . 
Notice that v = G Z£ for any (v; r) G Ti, that v G Zg has one or two nonzero entries for 
any (v; r) G T 2 , and that v G Z 2 J has zero, one, or two nonzero entries for any (v; r) G T 3 . 
Finally, let U be the set of all signed identities in Z 2 I S n . Notice that U C Z 2 I K. Recall 
that for any (u; e) G U, u G has exactly one nonzero entry. 

As with the classical model, before the two positions p and q have been chosen, the balls 
in the first rack may be permuted among themselves and the balls in the second rack may be 
permuted among themselves, without changing the state of the system following the switch. 
Thus, at each step, a random element of K is effectively generated whenever the procedure 
described above results in the identity or in a signed transposition in 7\; this event occurs 
with probability 

1 \n(n — 1) - r(n — r) n(n + 1) - 2r(n - r) 
2n 2n 2 4n 2 

A similar analysis shows that the procedure effectively generates a random element of 
UK = { (u; k) G Z 2 I S n : (u; e) G U and (0; k) G K} with probability a random element 
of T 2 K = {(v; tk) G Z 2 I S n : (v; r) G T 2 and (0; «) G K} with probability 3n(n-i)-6r(n-r) ^ 
and a random element of T 3 K = {(v;tk) G Z 2 I S n : (v;t) G T 3 and (0;k) G K} with 
probability 2r ^ a ~''} , 
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Notice that each element of T 3 if can be uniquely written as (v;tk), where (y;r) G T 3 
and (0;/c) G if. However, UK C T 2 if, with the exception that the elements (w; k) e UK 
with -u = (1, 0, . . . , 0) are not included in T 2 if when r — 1. But each element of UK can be 
uniquely written as (u; k), where (u; e) <E U and (0; k) G K and each element of T 2 K \ UK 
(where the set difference here is proper unless r = 1) can be uniquely written as (v; rr" 1 /?), 
where (v; r) G T 2 and (0; r _1 /t) G if. Let U r consist of the signed identities of the first r 
indices and U n - r consist of the signed identities of the last n — r indices. Notice that 

| If | = r\(n — r)\, 

\U r K\ = r r\(n — r)\, 

\U n - r K\ — (n — r)r\(n — r)\, 

\T 2 K \ UK\ = [\n(n - 1) - r(n - r)] r\(n - r)\, and 

\T 3 K\ = Ar(n — r) • r\{n — r)\. 

The signed Bernoulli-Laplace diffusion model may be modeled formally by a probabil- 
ity measure P on the hyperoctahedral group Z 2 I S n . We may thus define the following 
probability measure on the set of all signed permutations of Z 2 I S n : 

where (u; k) G U r K, 
where (u; k) G U n ^ r K, 



P(u; k 
P(v; tt~ x k 

P(v] TK 

P(x;ir 



2n 2 r\(n- 


-r)\ 


n—r 




2n 2 r\(n- 
1 


-r)\ 


2n 2 r\(n- 
1 


-r)! 


2n 2 r\(n- 


-r)\ 








(3.1.1) 

where (v; tt^k) G T 2 if \ UK, 
where (v; tk) G T 3 if , and 
otherwise. 



Since there are 2 n ■ n\ elements in Z 2 I S n , the uniform probability measure on the set of all 
signed permutations is given by 

U(x;tc) := for every (f; tt) G Z 2 I S n . (3.1.2) 

The following result establishes an upper bound on both the total variation distance and 
the £ 2 distance between P* k and U, where P* k is the probability measure on the homogeneous 
space X = (Z 2 I S n ) / (S r x S n - r ) induced by the convolution P* k of P with itself k times, 
and U is the uniform probability measure on X. It should be noted that \\P* k — C/||tv is the 
total variation distance to uniformity after k steps for the Markov chain on X defined above, 
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started at the chosen initial configuration. We establish an analogous result for the paired 
flips model as Theorem |3.4.3| . 



Theorem 3.1.3. Let P and U be the probability measures on the hyperoctahedral group 
Z 2 I S n defined in ( \3.1.J\ ) and ([3.1.Q ), respectively. Let P* k be the probability measure on the 



homogeneous space X = (Z 2 I S n ) / (S r x S n - r ) induced by P* k and let U be the uniform 
probability measure defined onX . Let k = |n(logn+c). Then there exists a universal constant 
b > such that 

||P**-I7||tv < \ [2 n - 0)] 1/2 ||^-^l|2 < be- c/2 for all c> 0. 



Notice that this is (essentially) the same result as that found in Theorem 2.5.3 . This Markov 
chain is twice as fast to converge as the random walk analyzed in Section 3 of Schoolfield 
(1999) in the special case of the hyperoctahedral group Z 2 I S n . 
In the following sections we present the results needed to prove this theorem and an 



analogous theorem for the paired flips model. In Section |3T2| we decompose the representation 
L(X) of Z 2 I S n into its irreducible constituents. In Section |3]3| we calculate the Fourier 
transform of the probability measure defined in ( |3.1.1| ), and this is followed by the proof of 



Theorem |3.1.3|. In Section FD| we perform a similar analysis of the paired flips model. 



3.2. Irreducible Representations in L(X). The decomposition of the representation 
L(X) of Z 2 I S n into its irreducible constituents is given by the following. See Section 3.4 
of Schoolfield (1999) for details about the irreducible representations of G I S n for any 67. 

Lemma 3.2.1. Let L(X) be the set of all (complex-valued) functionals defined on X = 
(Z 2 I S n )/(S r x S n - r ). Then, as a representation of Z 2 I S n , 

r n-r iAj (r—i)A((n—r)—j) 

L(X) = 000 P([(i+i)-<,fl;[n-(i+J)-m,m]) 

i=0 j=0 e=o m=0 

where pq( i+ j)^£ y e];[n-(i+j)-m,m}) i s the irreducible representation o/Z 2 I S n corresponding to the 
two-part partitions [(i+j) — £, i] and [n — {i+j) — m, m] ofi+j and n — (i + j) , respectively. 

Proof. Recall that as a representation of Z 2 I S n , L(X) is isomorphic to p Ts^s ™ ' 
where po is the trivial representation of S r x S n - r . It follows from Theorem 10 in Section 

*7 S S" S 

3.2 of Serre (1977) that p = p r ® Po"~ r where p r and p n ~ r are the trivial representations 
of S r and S n - r , respectively. 

Due to the transitivity of induction, 
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n t Z 2 lS„ _ f t (Z 2 ;SV) x (Z 2 ! s„_ r )l t z 2 ;5„ 

W I S r X S n - r - \P0 I 5 r X S n - r J I (Z 2 i S r ) X (Z 2 ! S n - r ) " 

Thus since {p* ®p 5 -} tg 2 x ( si* ^ ' ^ = R tg ? *} ® t£_ * 

follows that 

S r +Z 2 ; S r Sn-r + <Z 2 i S n -r I ■f-^lS 



M^J - \Po ls r ^ Po ls n _ r / l(Z 2 ; S r ) x (Z 2 ! S n _ r ) • 
It is a consequence of Corollary 4.4.7 of Greenhalgh (1989) that 

r n—r 

Po r \% 1 Sr = ©P(M;[r-*1) and Po" _r t£_'/"" r = ©P(b1;[(n-0-j]), 

i=0 j=0 

where P([i] ; [r— «]) is the irreducible representation of Z 2 ? SV corresponding to the trivial 
partitions [i] and [r — i] of i and r — i, respectively, and P([j] ; [(n-r)-j]) is the irreducible 
representation of Z 2 ? S^-r corresponding to the trivial partitions [j] and [n — r — j] of j and 
(n — r) — j, respectively. 

These results combine to show that 



s r n—r n 

L ( X ) ~ \ { P(mr-i])} ® {0P([7]i[(n-r)-j])} [ Tj 
^ i=0 j=0 ^ 



2 ? Sn 

2 2 i S r ) X (Z 2 I S„- 



r n—r 



= 9 9 {P(mr-i}) ® Pmi(n-r)-j])} tfzVi^) x (Z 2 I S n _ r ) ■ 

1=0 jr=0 

It follows from (the proof of) Lemma 4.4.5 of Greenhalgh (1989), which is a consequence 
of the "inducing- up rule" of Tokuyama (1984), that 

{P([i];[r-i]) ® P([M(n-r)-j])} ft) x (Z 2 )5„_ r ) P(W;M)> 

[A] M 

where the range of summation over partitions [A] of (i + j) is the range in the right-hand 
side of 

{PH ®Pbl} fsTx S, = 0PW 

[A] 

and, similarly, the range of summation over partitions \p] of n — (i + j) is the range in the 
right-hand side of 

{p[r-i] ®P[( n -r)-j}} fsZl^ S (n _ r) - 3 = P\A " 

M 

It follows from Corollary 4.4.7 of Greenhalgh (1989) that 
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{pih®pm} ts!" 



x S,- 







and 



{P[r-i\ ®P[{ n -r)-j]} T 



x 5, 



(n-r)— j 



(r-i)A((n-r)-j) 

P[n— (i+j)— m,m]- 

m=0 



These results combine to show that 

r n—r iAj (r—i)A((n—r)—j) 

L(X) = 000 P([(i+ j )-£ / ];[„-(i +i )- m , m ]). □ 



i=o j=o e=o 



m=0 



In Sections p.3| and |3.4j , it will be more convenient to use the following decomposition 
of L(X), which is a direct consequence of combining % and j and changing the order of 



summation in Lemma 3.2.1 



Corollary 3.2.2. Let L(X) be the set of all (complex-valued) junctionals defined on 
X = (Z 2 l S n )/(S r x S n - r ). Then, as a representation of Z 2 I S n , 

n Li/2J rA(j-l) (r-i)A((n-j)-(r-i)) 

L(X) = 00 P{]3-lA^3)-mM) 

j=0 e=o i=iv(r-(n-j)) m=0 

where p{\j^^.,[{n-j)-m,m}) ^ s the irreducible representation of Z 2 I S n corresponding to the 
partitions [j — £, £} and [n — j — m, m] of j and in — j), respectively. 

The number of times that a particular representation p occurs in the direct sums in 
Lemma 3.2.1 and Corollary 3.2.2| is its multiplicity m p in the decomposition of L(X). 



3.3. Analysis of the Independent Flips Model In order to apply the Upper Bound 
Lemma Q2.2.4 ), we must now calculate the Fourier transform at each nontrivial irreducible 
representation of Z 2 I S n occurring in the decomposition of L(X). We use the same technique 
as was used in Section 1273. 



Lemma 3.3.1. Let P be the probability measure on Z 2 I S n defined in ( \3.1.1\ ). Let 
P = /°(b'-^];[(n-i)-m,m]) be the nontrivial irreducible representation o/Z 2 I S n , corresponding 
to the partitions [j — £, £] and [(n — j) — m, m] of j and n — j , respectively, and occurring in 
the decomposition of L(X). Then, in a certain basis, the Fourier transform is 

2£(j-£+l) 



P{P) 



J 



rr 



I(P) 
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where I(p) is the d p x d p matrix 



10 
00 



and I is the m„- dimensional identity matrix. 



Proof. Recall that the trivial representation of Z 2 I S n corresponds to the partition 
[n] of n. Thus we must calculate the Fourier transform for the other nontrivial irreducible 
representations of Z 2 I S n occurring in the decomposition of L(X), which were found in 
Corollary ggg . Notice that 

p( P )= p@;k)p&,k) + J2 + J2 

(5;«)eK 



(u;K)£U r K 



(u;K,)£Un- r K 



+ P(v;tt 1 k)p(v;tt x k) + P(v; tk)p(v; tk). 

{v;tt- 1 k)£T 2 K\UK (v;tk)€T 3 K 

Choose an ortho normal basis in V such that the first m p basis vectors are K- invariant, as 
described in Section |2.2| . It then follows from Lemma |2.2.3| that, in this basis, 



E,- n ,- s wfn + 1) — 2r(n — r) v—v 



(0;k)gK 



(0;«)eJC 



n(n + 1) — 2r(n — r) 
An 2 

n(n + 1) — 2r(n — r) 
An 2 



r\(n — r)\ 



(0;k)6A' 



Recall that C/X C T 2 ^- This is a result of the fact that for two-thirds of the elements 
in T 2 , namely the odd signed transpositions, v G Zg ^ as exactly one nonzero element. Let 
(ui] e) G U be the signed identity whose only nonzero entry in occurs in position i. Let 
(vi] T{j) G T 2 be the odd transposition which transposes the entries in positions i and j and 
whose only nonzero entry in Vi occurs in position i. 

Notice that for any k G K, e)(0; k) = (uf, k) = (vf, t^t^k) = (vl; Ty)(0; k) for 
i 7^ j. So it follows that 

(u;K,)£U r K 



2n 2 r\(n — r)\ 



(u;K)eU r K 



2n 2 r\(n — r)\ ' 
v ' i=l 



E E 



(0;k)gA- 



1 r 

-T 



2n 2 



8=1 



r\(n — r)\ 



(0;/t)eA- 
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Similarly, 



Y p & K )p(u; k) 



1 - 

2tf S 



i=r+l 



'J / 



r!(n — r)! 



(o ; «;)ei<r 



Thus, since 



P(v;tt 1 k)p(v;tt 1 k) 

{v;tt- 1 k)&T 2 K\UK 



1 

2^ 



run — r 



.(t?;r)eT+ 

where T 2 + is the set of all even transpositions in T 2 , it follows that 
E P(u;k)p(u;k) + Y P (u; k) p(u; k) 

(u;K.)£U r K (u\K)(^U n - r K 

+ Y^ P(tf; rr _1 fi;)p(-u; tt^k) 

(v;tt- 1 k)€T 2 K\UK 



Y p& k ) 



1 

2^? 



{u;e)€U (v;t)£T 2 



r\(n — r)! 



Furthermore, we have 

E P(v;tk)p(v;tk) = 

(v;tk)£T 3 K 

These results combine to show that 



2n 2 



P{p) 



n{n + 1) — 2r(n — r) 



An 2 



Y p(^ t ) 



Hp) 



r\(n — r)\ 



Y 

(0;k)£K 



Y 



1 

2^2 



E 

(u;e)e!7 



E 



p(u;e) + V p(V;r) - 



(tf;r)eTi 

It was shown in Section 3.3 of Schoolfield (1999) that T splits into two conjugacy classes 
in Z 2 I S n , namely, the even transpositions (v; r + ) (which change neither or both of the 
charges) and the odd transpositions (v;t~) (which change exactly one of the charges). It 



then follows from Lemma 2.4.1 that 
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£ pftr) = [§|T| r+(p) + \\T\r-{p)] I = n(n - 1) [r+(p) + r~(p)] I, 

where r + (p) := Xp(v',T + )/d p , r~(p) := xp(v'i T ~)/d P , and / is the c^-dimensional identity 
matrix. It follows from Lemma 3.5.1 of Schoolfield (1999) that 

j(j-l)r(X) + (n- j)(n- j -l)r(fi) 



r + ( P ) 



r (p) 



n{n — 1) 

j(j-l)r(X) - (n- j)(n-j - l)r(fi) 



n(n — 1) 

where r(A) := X[A](r)/d[A], with [A] = [j - £,£] and r e and r(p) := XmCOM,*], witn 
[fit] = [n — j — m, m] and r e SVi-j- From these results it follows that 

£ P (t?;t) = 2j(J - 1) r(A) I. 

(v;t)£T 

Since = po Tx ' ™> where po is the trivial representation of K, it follows that for 

any p occurring in the decomposition of L(X), p J.^ 2 1 Sn is the direct sum of d p copies of po- 
Thus since T\ C if, we have 



[|n(n — 1) — r(n — r)] /. 



(■u;t)€Ti 

It was shown in Section 3.3 of Schoolfield (1999) that U is a conjugacy class in Z 2 I S n . 
It follows from Lemma |2.4.1| that 

£ P(^ e ) = 1^1 [Xp(u;e)/d p } I, 

(n;e)e(7 

where the right hand side is calculated at any (u; e) G U. It follows from Lemma 3.5.1 of 
Schoolfield (1999) that Xpi^'i e )/d p = (2j — n)/n. From these results it follows that 



(u;e)eU 

These results combine to show that 

n(n + 1) — 2r(n — r) 



(2j - n) I. 



P(P) 



An 2 



Hp) 



i 



+ — [(2j - n) + 2j(j - 1) r(A) - §n(n - 1) + r(n - r)] J(p) 



2n 2 

J i(i - 1) 



^ + 



r(A) 
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Recall that it follows from Lemma 12.4.21 that 



r(A) 



l(j-t)(j-£-l) + i(i-S)] 



Therefore, for the nontrivial irreducible representation of Z2 I S n , corresponding to the 
partitions [j — £, £] and [(n — j) — m,m] of j and n — j, respectively, and occurring in the 
decomposition of L(X), 



P{p) 



j +U-£)U-£-i) + t(t-$ 



ir 

_ 2£(j-e + i) 

n 2 ri 1 



Hp) 



Hp)- □ 



We have now established the results necessary to prove Theorem 3.1.3 



Proof of Theorem |3.1.3| . Notice that the probability measure P defined in Q3.1.1 
is clearly right i^-invariant. (In fact, it is bi-i^-invariant.) By applying the result from 
Lemma |3.3.1| to the Upper Bound Lemma ( 2.2.4|) , we find that 



\P* k -U\\ 2 TY < \ [2 n ■ (™)]||P* fc - U\ 



P m p 



f 2 £(j-£ + l) 



-i 2k 



rr 



where the sum is taken over all nontrivial irreducible representations p = P([j~£/]-[( n -j)-m,m]) 
of Z2 I S n occurring at least once in the decomposition of L(X). The factor m p comes from 
the trace of I(p). It follows from Lemma 3.5.1 of Schoolfield (1999) that d p = (")<iy_^^] ■ 
d[(n-j)- m ,m]- ^ ^en follows from Lemma [2.3. 1| that 



d P = [© - G-i)] ■ - cri)] 



For notational purposes, define 



F{j,- 



2£{j-£+l) 



Rather than explicitly calculate m p , we will sum over the entire range of indices in the 
decomposition of L(X) found in Corollary [3.2.2 , thereby including each representation p a 
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total of m p times in the summation. Thus we have 
\\P^-Uf TV < |[2"-C)]||P^-^ 

n U/ 2 J rA(j-£) (r-i)f\((n-j)~(r-i)) 

= iEE E E G) [0 - (A)] ■ [( V) - Cri)] • *u *) 

j=0 £=0 i=£V{r-(n-j)) m=0 



2k 



Li/2J 

iEC)E r// 

Li/2J 

iE©E r ' 

rA(j-i) 



r/\(j-l) (r-i)A((n-j)-(r-i)) 

2k E E K^)-(::i)] 

i=£V(r-(n-j)) m=0 



rA(j-e) 
i=£V(r-(n-j)) 



(3.3.2) 



Notice that C-f) ^ E (V) = 2 "~ J - Thus we have 

i=£V(r-(n-j)) «=0 



Li/2J r 



<iE 2 

3=0 



n-j 



E 

£=0 



J 



f 2£(j-£ + l) 

2 



-i 2k 



rr 



For each 1 < j < n, it follows from (the calculations in the proof of) Theorem 2.5.3 that 
we may bound all but the £ = term of the inner sum above by 

Li/2J r / -\ / \ ~i r -9 nnl ■ „ , ,M2k 



E 



J 



1 



f 2£(j-£+l) 

2 



n 



4fc Li/2J 

n I E/ 



2£(j-£ + l) 



n 2fc 



-Ik 



< ( - I 4a z e 
n 



2„-c 



for a universal constant a > 0, when /c > |j(log j + c). Since n > j, this is also true when 
k > ^n(logn + c). 

We must also bound the term for the trivial representation [A] = [j] for 1 < j < n — 1. 
Since in these cases £ = 0, we have 



J 

- 1 



f 2£(j-£ + l) 

2 



n 



-i 2fc 













These results lead to the upper bound 
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\\p* k -ur TV <i[2-.(-)]\\p^-u\\ 



3=0 



3=1 



n 
3 

n-j 



IJ/2J 

E 

J/ V™ 



J 2 2£(j-£ + l) 



n- 



2A- 



n-l 



J J \n 



Ik 



Now notice that, when k = \n(\ogn + c), then 



4[3 nl °g( n )+3 cn ] f j\ -n[-log(n)-c] 



?7 



which combines with the results above to give 

||? :fe -^|| 2 :v<i[2«-P]||P^-/7|| 



<a 2 e" c ^2 

n-l 



n-j 



n 



3 J \ n 



+ iE 2 ^' 

If we let i = n — j it follows that 

||P* - m v < \ [2" ■ (»)] ||P*-tf||: 



7/ 



3 J \ n 



-nlogl 



-nlogl £ . 



n-l 



i=0 
n— 1 



n-l 



i=i 

n-2 



a 2 e" c 



(z + l)! 



< a 2 e" c exp (2e" c ) + §e~ c exp (2e~ c ) . 
Since c > 0, we have exp(2e~ c ) < e 2 . Therefore 

||P^-^|| 2 V < i[2«-C')]||P^-t/|| 2 < [(a 2 + |)e 2 ]e" c , 

from which the desired result follows. □ 

Theorem 3.1.3| shows that k = \n (log n + c) steps are sufficient for the (normalized) 
£ 2 distance, and hence the total variation distance, to become small. A lower bound in 
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the (normalized) £ 2 metric can also be derived by examining 2n (l — -^) 4k , which is the 
dominant contribution to the summation ( 3.3.2 ) from the proof of Theorem [3.1.3| . This term 



corresponds to the choice j 



n 



1 and t = 0. Notice that k = \n (log n — c) steps are 



necessary for just this term to become small. 

Recall that Theorem 2.5.5 shows that, for values of r not too small, k — \n (logn — c) steps 
are necessary for the total variation distance to uniformity to become small in our variant 
of the Bernoulli-Laplace diffusion model. This is exactly the independent flips model, if the 
charges of the balls are ignored. For such r (in particular, for the symmetric case r = n/2), 



Theorem 2.5.5 provides a matching lower bound on the distance to uniformity in the total 



variation metric. The upper bound in Theorem |3.1.3| , taken together with its matching lower 
bound, gives another example of the "cutoff phenomenon." 



3.4. Analysis of the Paired Flips Model. At each step of the process introduced in 
Section |3J], the charges of the switched balls were changed independently. Suppose instead 
that, at each step, the charges of the switched balls are either both changed or both not 
changed. It is this type of process that we now examine. 

We now describe a slight variant of the independent flips model introduced in Section |3.1| . 
Suppose that the balls and racks are as described in Section |3J]. At each step, independently 
choose two positions p and q uniformly from {1,2, ... ,n}. 

If p ^ q, switch the balls in positions p and q. Then independently, with probability 
1/2, change the charge of the ball moved to position p and change the charge of the ball 
moved to position q. Notice that this combination of operations is an even transposition; the 
probability that an odd transposition occurs now vanishes. Then, if necessary, permute the 
balls on each of the two racks so that their labels are in increasing order. 

If p = q (which occurs with probability 1/n), leave the balls in their current positions. 
Then, again independently with probability 1/2, change the charge of the ball in position 
p = q. The probabilities of the identity and of the signed identities are thus unchanged 
from the independent flips model. We refer to the process described above as the paired flips 
model. 

The following analysis, while similar in format, is different in content from that found in 



Section |3.1| . Let T be the set of all even transpositions in Z 2 I S n . Also, let T\ be the set of 
all even transpositions in K, let T 2 be the set of all even transpositions in (Z 2 I K) \ K, and 
let T 3 be the set of all even transpositions in (Z 2 I S n ) \ (Z 2 I K). Thus T = 7\ U T 2 U T 3 . 
Notice that v = 6 Z 2 for any (v; r) E Ti, that v E Z 2 has exactly two nonzero entries for 
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any (v;t) E T 2 , and that v E Z£ has zero or two nonzero entries for any (v;r) E T 3 . Finally, 
let U be the set of all signed identities in Z 2 I S n . Notice that U C Z 2 I K. Recall that for 
any (u; e) E U, u E Z 2 has exactly one nonzero entry. 

As with the classical model, before the two positions p and q have been chosen, the balls 
in the first rack may be permuted among themselves and the balls in the second rack may be 
permuted among themselves, without changing the state of the system following the switch. 
Thus, at each step, a random element of K is effectively generated whenever the procedure 
described above results in the identity or in an even transposition in Ti; this event occurs 
with probability 

1 \n(n — 1) — r(n — r) n 2 — 2r(n — r) 
2n n 2 2n 2 

A similar analysis shows that the procedure effectively generates a random element of UK 
with probability a random element of T 2 K with probability ra(ra ~ 1 - > 2 ~^ r(w ~ r - ) , and a random 
element of T%K with probability 2r "( n ~ r ) . 

Notice that each element of T%K can be uniquely written as (v; tk,), where (v; r) E T3 and 
(0; k) E K, each element of T 2 K can be uniquely written as (v; rr _1 /t), where (v; r) E T 2 and 
(0; t~ 1 k) E K, and each element of UK can be uniquely written as (u; n), where (m; e) E U 
and (0; k) E K. Thus 

I If I = r!(n — r)!, 
\UK\ = nr\(n — r)!, 

\T 2 K I = [|n(n — 1) — r{n — r)] r!(n — r)!, and 

|T 3 if| =2r(n-r) -r!(n-r)!. 

This paired flip variant of the signed Bernoulli-Laplace diffusion model may be modeled 
formally by a probability measure Q on the hyperoctahedral group Z 2 I S n . We may thus 
define the following probability measure on the set of all signed permutations of Z 2 I S n : 

3(0; «) := fegg^j where (0; «) G K, 

Q($> K ) 2n^r!(n-r)! Wllere K ) ^ U K , 

Q(v; rr _1 K) := ra2r | ( i t _ r) | where (u; rr" 1 ^ G T 2 if , (3.4.1) 
tk) := Tt 2 r!(n— r)! where ({7; tk) E T 3 K, and 
(^(a 7 ; 7r) := otherwise. 
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In order to continue our analysis of the paired flips model, we must now calculate the 
Fourier transform of Q at each nontrivial irreducible representation of Z 2 I S n occurring in 



the decomposition of L(X). We use the same technique as was used in Section [O 



Lemma 3.4.2. Let Q be the probability measure on Z 2 I S n defined in ($-4- A) - Let 
P = P{[j-^A;[(n-j)~m,m]) be the nontrivial irreducible representation o/Z 2 I S n , corresponding 
the partitions [j — £, £] and [{n — j) — m, m] of j and n — j , respectively, and occurring in the 
decomposition of L(X). Then, in a certain basis, the Fourier transform is 



Q{p) 



21Q-1+1) , (n-j) 2 
n 2 7)2 



2m(n—j—m+l) 

7a 



n-3 

,„2 



Hp) 



where I(p) is the d p x d p matrix 



10 
00 



and I is the m p - dimensional identity matrix. 



Proof. Recall that the trivial representation of Z 2 I S n corresponds to the partition 
[n] of n. Thus we must calculate the Fourier transform for the other nontrivial irreducible 
representations of Z 2 I S n occurring in the decomposition of L(X), which were found in 
Corollary gTgTgj . Notice that 



Q{p)= Q&k)p(0;k) + Yl Q(u;k) P (u;k) 

+ ^2 Q(v;tt~ 1 k)p(v;tt~ 1 k) + Q(v;tk)p(v;tk). 

{v;tt- 1 k)&T 2 K {v;tk)&T 3 K 

Choose an ortho normal basis in V such that the first m p basis vectors are .fT- invariant, as 



described in Section |27|. It then follows from Lemma |2.2.3| that, in this basis, 

n 2 — 2r(n — r) 



Q(0;k) P (0; 

(0;k)G/< 



2n 2 r\(n — r) 

n 2 — 2r(n — r) 
2^2 

n 2 — 2r(n — r) 
2m 2 



E p(0;«) 



r!(n — r)\ 

Tip)- 



Notice that 
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E K M 



u: k) 



{u;k)€UK 



E 

(v;tt- 1 k)£_T 2 K 



1 

2^ 



E * 

(u;e)GC/ 



u: e 



r\(n — r)! 



E 



Q(v;tt 1 k)p(v;tt 1 k) = — 



rr 



E o< 

(u;-r)eT 2 



r!(n — r)! 



E ^ 



and Q(v]tk)p(v]tk) = — 

{v;tk)€T 3 K 

These results combine to show that 
n 2 — 2r(n — r) 



E P^]t) 

(v;r)eT 3 



r\(n — r)\ 



E p& 



K 



Q(p) 



2n 2 



Hp) 



i 

+ — 

vr 



\ E p& e ) + E p^i t ) - E ^ ;r ) 

(u;e)eE/ (v;r)eT (^r)eTi 



Recall that the even transpositions T form a conjugacy class in Z 2 I S n . It then follows 



from Lemma 2.4.1 that 



E p(^ t ) 



= \T\ r (p) I = n(n — 1) r (p) I, 
where r + (p) := Xp(v'i T+ )/d p for any even transposition (v; r + ) and J is the ^-dimensional 



identity matrix. Recall that r + (p) was determined in the proof of Lemma |3.3.1| . From these 
results it follows that 

E P&t) = [j(j-l)r(X) + (n- j)(n-j- 1) r{p)\ /, 

(v;t)&T 

where r(A) := x [A ](r)/rf [A] , with [A] = [j - £,£) and r G S,-, and r(//) := Xm( t )Mm] with 
[//] = [(n — j) — m, m] and r G 5* n _j. As in the proof of Lemma [3.3. 1| , 

p(v;r) = \Ti\ I = [|n(n — 1) — r(n — r)] / and 

(o;t)6Ti 

E ^("5 e ) = _ n ) L 

These results combine to show that 
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Q(p) 



n — 2r{n — r) 
2^ 



I(p) 



+ 



^(2j-n) + j(j-l) r(A) + (n-j)(n-j-l) r(n) - \n{n-l) + r(n-r) 



Hp) 



n 



!) /u (n-j)(n-j - 1) , ' 
2 — r (A) + 2 r/i 



Recall that it follows from Lemma |2.4.2| that 

r(A) = — L_r(.? - ^) 1) + ^-3)] and 



r(//) 



1 



r \( n — j — rn)(n — j — m — 1) + m(m — 3)1 . 

(n-j)(n-j-l) 

Therefore, for the nontrivial irreducible representation of Z 2 I S n , corresponding the par- 
titions [j — £,£] and [(n — j) — m,m] of j and n — j, respectively, and occurring in the 
decomposition of L(X), 



Q{p) 



j +0'-^)0'-^-l) + 1(1-3) + (n-j-m)(n-j-m-l) + m(m-3) 

n.2 



ggj^W) , (n-j) 2 



2m(n- j— m+1) 



„2 



Hp) 
Hp)- 



□ 



The following result establishes an upper bound on both the total variation distance and 
the £ 2 distance between Q* k and U, where Q* k is the probability measure on the homogeneous 
space X = (Z 2 I S n )/(S r x 5 n _ r ), induced by the convolution Q* k of Q with itself k times, 
and U is the uniform probability measure on X. 

THEOREM 3.4.3. Let Q and U be the probability measures on the hyperoctahedral group 
Z 2 I S n defined in ( $.4-l\ ) and ( 3.1.2J , respectively. Let Q* k be the probability measure on the 
homogeneous space X = (Z 2 I S n ) / (S r x S n - r ) induced by Q* k and let U be the uniform 
probability measure defined onX . Let k = |n(logri + c). Then there exists a universal constant 
b > such that 

WQ^-Uhv < l[2".^)] 1/2 \\qr k -U\\ 2 < be- c l 2 for all c> 0. 
Proof. Notice that the probability measure Q defined in ( |3.4.1|) is clearly right K- 



invariant. (In fact, it is bi-K- invariant.) By applying the result from Lemma 3.4.2 to the 



Upper Bound Lemma ( 2.2.4 ), we find that 
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\\Q* k -U\\ 2 TY < J[2*.Q]||Q**-[/|| 

= i d p m t 



21Q-1+1) , (n-j) 2 

n.2 n.2 



2m(n-j-m+l) 

r) .2 



n 2fc 



r>2 



where the sum is taken over all nontrivial irreducible representations p = P([j-^];[(n-j)-m,m]) 
of Z2 ? occurring at least once in the decomposition of L(X). Recall from the proof of 



Theorem 13.1.31 that 



(■)[©- (/,)]■[»-(::!)]■ 



As in the proof of Theorem |3.1.3| , rather than explicitly calculate m p , we will sum over 
the entire range of indices in the decomposition of L(X) found in Corollary |3.2.2| , thereby 
including each representation p a total of m p times in the summation. Thus we have 



\\Q* k -u\\l w < l[ 2n -C)]\\Q* k - u \\\ 

E © [ffl - ■ kv) - (.r-i)] (3.4.4) 



n Li/2J rA(i-f) (r-i)A((n-j)-(r-i)) 

IE E E 

j=0 £=0 i=A/(r-(n-j)) m=0 



2£(j-£+l) , (n-j) 2 
9 ~~ r~ 9 



Notice that, when 1 < j ' < n — 1, 



2l(j-l+l) , (n-j) 2 
9 "t - 9 



2m(n—j — ra+ 1) 



2m(n— m+1) 
^2 



2fc 



r,2 



2 A- 



2m(n— j— m+1) 1 
(n-j) 2 n-j 



< max J (J) 



< ta 



- 2fe 



2l(j-l+l) 



2A- 



2£(j-e+i) 



3 

2k 



2k 



2 k 



1 - 



2m(n— j — m+1) 



, 2k 



1 - 



2m(n— j— m+1) 1 



2fc 

1 

n-j 
2k 



2k 



(n-j) 2 



where the first inequality is due to the fact that (ax + (1 — a)y) 2k < max jx 2fc , y 2fc }, for 
< a < 1. It can be determined that, for all possible choices of j and m with 1 < j < n — 1 
and m < (n — j)/2, 

j _ 2m(n-j-m+l) _ 1 ^ q 



(n-j) 



n-j 



except when j = n — 2 and m = 1. But when j = n — 2, notice that 



36 



C.H. SCHOOLFIELD, JR. 



2m(n—j—m+l) 1 
(n-j) 2 n-j 



2k 



< 



E[ 



: 2m(n— j— m+1) 

1 FRF 



m=0 m=0 

These results combine to give the upper bound 



2fc 



IIQ* fc -^llTV < J I 2 " •OHIO**- ^11 



^ 14 



n [j/2\ rA(j-e) (r-i)A((n-j)-(r-i)) 

i£E E E (")[©-(/■)]■[«- (::!)] 

j=l £=0 i=A/(r-(n-j)) m=0 



X 



(i) 



2/,' 



1 - 



2i(j-e+i) 



21,' 



n-1 [j/2\ rA(j-t) (r-i)A((n-j)-(r-i)) 

+ \i E E E E © [(8 - ((-J] ■ [(",;') - (::'■)] 

jr=0 £=0 j=A/(r-(n-j)) m=0 



X 



2A; 



1 - 



2m(n— j — m+1) 

0-j) 2 



21: 



, 1 / n-l \ 2fc 
4 V n / ' 



where we must modify the first iterated summation to exclude the term for j = n and i = 
(and hence i = r and m — 0) and we must also modify the second iterated summation to 
exclude the term for j = and m = (and hence £ = and « = 0). The final expression 
reintroduces the appropriate term in the exact formula for the squared £ 2 distance at the 
representation corresponding to j = and m — 0. 

Notice that if, in the second set of braces | | , we put i' — r — i and j' = (n — r) — j and 
change the order of summation, we obtain 



IIQ* fc -^llTV < i[2 n -C)]ll^ fc -^ll 



< 



n [j/2\ rA(j-i) (r-i)A((n-j)-(r-i)) 

i E E E E K9 - (A)] • [("„?) - (r-i)l 

J=l £=0 i=N(r-(n-j)) m=0 



X 



(i) 



2A- 



i"2 



2/,: 



n-1 L(™-i)/ 2 J rA((n-j)-m) (r-i)A(j-(r-i)) 

+ UE E E E [© - (A)] • KV) - (n)] 

j=0 m=0 i=mV(r-j) ^=0 



X 



2A; 



2m(n— j— m+1) 
1 RTF 



2 k 
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Notice that if we now put j' = n — j and interchange the roles of £ and m, then the second 
summation becomes identical to the first. Thus, combining these summations and continuing 
as in the proof of Theorem |3.1.3| , we have 



\\Q* k -u\\ 2 TV < \[r>C)]\\Q*-u\\ 



< 



n [j/2\ rA(j-i) (r-i)A((n-j)-(r-i)) 

IE E E E (")[©- • K"^) - (;:!)] 

j=l £=o i=(.y(r-(n-j)) m=0 



X 



(a 



2k 



2e(j-£+l) 



2k 



4 V n J 



1 \2k 



< 5E 2 ' 

3=1 



n-:i n \ 3 



2k U/2J 



2£(j-£+l) 



2/,- 



+ 



n — 1 



n 



2A- 



Recall from the proof of Theorem 3.1.3| that, when k > \n(\ogn + c), we may bound the 
inner sum above using 



2k Li/2J 

E 



for 1 < j < n — 1 , and using 

2k U/2\ 

E 



1 - 



J 

- 1 



2£(j-£ + l) 



~1 2fc 



,;2 



< 



2£(j-£+l) 



~1 2k 



2k 



Aa 2 e~ c + 1 3 
n 



< 4aV c 



2A- 



for j = n. So this is also true when k = |n(logn + c). 
Now notice that when k = |n(logn + c) 



\ 2k /„-c\ -nlogO'/n) 
J \ I c 



n 



n 



These results lead to the upper bound 

\\Q*-u\\h<\[^'C)] \\Q*-v\ 
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Continuing as in the proof of Theorem |3.1.3| , we find that 
W k -U\\\ Y <\[2--{-)] \\Q^-U\\l 

< 2aV c exp (2e- c ) + e~ c exp (2e~ c ) + \e~ c 
Since c > 0, we have exp(2e _c ) < e 2 . Therefore 

IIQ^-^IItv < W n -0]W k -U\\l < [(2a 2 + l)e 2 + i]e- 
from which the desired result follows. □ 



Theorem |3.4.3| shows that k — |n (log n + c) steps are sufficient for the (normalized) 
£ 2 distance, and hence the total variation distance, to become small. A lower bound in 
the (normalized) £ 2 metric can also be derived by examining 2n (l — ^) 4A \ which is the 
dominant contribution to the summation (|3.4.4j ) from the proof of Theorem |3.4.3| . This term 



corresponds to the choice j = n — 1 and £ = 0. Notice that k = \n{\ogn — c) steps are 
necessary for just this term to become small. Furthermore, our upper (|nlogn) and lower 
(^n \ogn) bounds on the number of steps required for the (normalized) £ 2 distance to become 
small differ by a constant factor. We have not been able to close this gap. 

Recall that Theorem [2.5.5| shows that, for values of r not too small, k — \n (logn — c) 



steps are necessary for the total variation distance to uniformity to become small in our 
variant of the Bernoulli-Laplace diffusion model. This is exactly the paired flips model, if 
the charges of the balls are ignored. For such r (in particular, for the symmetric case r = n/2), 
Theorem |2.5.5| provides a lower bound (differing from the upper bound by only a constant 



factor) on the distance to uniformity in the total variation metric, just as in Section |0 
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